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Auroral Micropulsation Instability 
F. V. Coroniti* and C. F. Kennel 

Physics Department, University of California, Los Angeles, California 90024 


Abstract 

In this paper we describe a drift instability of Alfven 
waves driven by the sharp electron thermal gradient at the 
inner edge of the electron plasma sheet. The analysis predicts 
a localized instability with wavelengths perpendicular to the 
magnetic field comparable with an ion cyclotron radius and 
parallel wavelengths comparable with the length of field lines. 
The fastest growing wave should have an - 10 second period, 
consistent with observations of Pi 1 micropulsations observed 
during substorms. 


* Also at Physics Department, and Space Science Laboratory, University 


of California, Berkeley, California. 



- 1 - 


1. Introduction 

Intense low frequency micropulsations occur during magnetic storms 
and auroral substorms. Ground magnetometer observations (McPherron, et 
al., 1968) reveal a variety of oscillations, each with its own locah time 
and spectral characteristics > whose complex morphology makes it unlikely 
that a single mechanism explain all the observations. We attempt to 
identify - 10 second period magnetic micropulsations, which are morning 
side manifestations of substorms associated with modulations of energetic 
electron precipitation fluxes (Parks, et al., 1968a) with a particular 
drift instability of Alfven waves, driven by the strong electron thermal 
gradient at the inner edge of the electron plasma sheet (Vasyliunas, 

1968a, b; Frank, 1968). 

Except for the Pc-1 band, micropulsation frequencies are much smaller 
than the equatorial proton cyclotron frequency. Thus, except for Pc-1, 
only extremely high energy particles can be in cyclotron resonance, of 
which there are presumably too few to trigger an instability. We therefore 
rule out cyclotron resonance interactions and pitch angle anisotropy as a 
low frequency micropulsation source. 

Multiplying the long wave periods by a characteristic propagation 
speed (such as the Alfven speed) suggests that at least one micropulsation 
spatial dimension should be roughly the length of the line of force. Yet 
the association of 5-40 second period micropulsations with precipitation 
pulsations (McPherron, et al., 1968) yields indirect evidence that 
micropulsations are localized transverse to the magnetic field. Due to 
our poor understanding of mode structure, propagation, and ionospheric 
coupling, ground magnetometer measurements do not fix reliably the size and 
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location of micropulsation generation regions. Since electrons are guided 
by the magnetic- field, electron precipitation pulsations give inherently 
better information than micropulsations, which cannot be guaranteed 
a priori to be so guided. Barcus and Rosenberg (1965), Barcus, et al., 
(1966), and Parks, et al., (1968a) , found electron precipitation pulsations 
to be limited to 100-200 km regions in the ionosphere, which project onto 
the equatorial plane as several thousand kilometer regions, less by far 
than the length of auroral lines of force. 

The large scale lengths parallel and small lengths transverse to the 
magnetic field direction, together with the low frequency, are consistent 
with drift waves. Furthermore, drift instabilities derive their energy 
from spatial gradients in the distribution of resonant particles. In 
straight magnetic fields, the important wave-particle interaction is Landau 
resonance (i.e., the particle parallel velocity v. matches the parallel 
phase velocity w/k^ ); the analog of Landau resonance in mirror geometry 
is bounce resonance. Cyclotron resonances do not affect low frequency 
drift instabilities. 

The steep inner edge of the electron plasma sheet (Vasyliunas, 1968a) 
is an obvious candidate for various drift instabilities. We search here 
for one ttfhich produces magnetic micropulsations between the ion and 
electron bounce frequencies. While magnetic fluctuations on the ground 
could originate from a rotating ionospheric line current driven by a 
purely electrostatic wave (Wilson, 1966), a simpler hypothesis is that 
micropulsations are magnetically polarized in deep space. We therefore 
generalize, for magnetospheric conditions, an analysis by Mikhailovsky 
and Rudakov (1963) of drift instabilities of Alfven waves. D'Angelo (1969) 
has also suggested drift instabilities are responsible for micropulsations. 
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Earlier, Swift (1967) analyzed an electrostatic drift instability without 
considering coupling to an oscillating magnetic field, 

Section 2 elucidates when Alfven waves can have drift couplings to 
resonant particles, using physical arguments drawn from the two- fluid 
theory of low frequency waves in spatially homogeneous plasmas and a 
uniform magnetic field. In section 3, we derive a kinetic theory dispersion 
relation for drift waves in a spatially inhomogeneous plasma, assuming the 
magnetic field approximately uniform. While the techniques used are 
standard, many approximations are needed. These are discussed in 3.2, 
and the calculation described in 3.3. Those uninterested in calculational 
details may turn to section 4, where a simple dispersion relation, 
appropriate to the auroral electron boundary, is described. Here the range 
of unstable frequencies and wavelengths for Alfven waves in a strong 
temperature gradient is outlined. A comparison with observations in 
section 5 indicates that the present ideas are consistent with the 
observed spatial localization, magnetic polarization, and range of unstable 
frequencies and wavelengths of 10 second micropulsations. However, further 
theoretical work removing many of our simplifying assumptions is needed 
for truly definitive conclusions. 
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2 . Motivation of Drift Approximation 
2.1) Introduction 

Information about changing plasma conditions propagates at various 

characteristic wave and particle speeds. For example, in single-fluid 

hydromagnetics, the relevant waves are the fast, intermediate, and slow 

waves (Kantrowitz and Petschek, 1966). In spatially inhomogeneous plasmas, 

particle drifts are additional physically significant propagation speeds. 

When the magnetic field strength B and the ion (+) or electron (-) 

+ 

pressure P” have spatial gradients, the corresponding drifts are 


+ vj B x VB 




B x V in P + 


P" B 


( 2 . 1 ) 


where v is the particle's velocity component perpendicular to the 
magnetic field, p* is the ion or electron mass density, B is the 
magnetic field, and = ± eB/M + c is the ion or electron frequency. 

( e is the magnitude of the electronic charge in esu; c , the velocity 

of light; and M + , the ion or electron mass. Henceforth, we will use 

+ + 

Gaussian units.) v Q is an actual guiding center drift, while v p , 

4 - -■ 

the pressure drift, is not; v D - v_ is the macroscopic drift between 

rt,* '•x.r 

electrons and ions providing the diamagnetic magnetic field gradients 

associated with pressure gradients. We eventually treat B to be 

+ 

approximately uniform; therefore we concentrate on v p . Let us now estimate 
its order of magnitude: 




+ 


f 


V 

p 


( 2 . 2 ) 
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where is the ion thermal speed, R + = c , the thermal ion cyclotron 

+ 4 ; 

radius, and T the ion or electron temperature in energy units. L is a 

4 4 

logarithmic pressure scale length. Ordinarily R + /L^ « 1 , so that v^~ is 

much smaller than the particle thermal speeds. 

For weak spatial inhomogeneities, the dispersion relation obtained 
from homogeneous plasma theory often describes wave properties well. One 
index of the importance of drift corrections to the dispersion relation is 

4 

the ratio of v^ _ , which is perpendicular to the magnetic field, to the 
perpendicular phase velocity obtained from homogeneous theory. When 
k^VpVco ~ 1 , plasma inhomogeneity significantly modifies the wave propagation 

4 

speeds and unstable growth rates. Since v^ is small, only small 
perpendicular phase velocity waves have important drift corrections. We 
search for such waves in section 2.2. 


2.2) Low Frequency Waves in the Two-Fluid Approximation. 

By treating the ion and electron gases separately, the two-fluid 
description allows for the effects of the finite inertia of the individual 
particle species, which manifests itself in the inability of first ions 
and then electrons to produce rapid variations in plasma current. The 
two-fluid dispersion relation for the three wave modes with frequencies 
below the electron plasma frequency is (Braginskii, 1957; Stringer, 1963; 
Formisano and Kennel, 1969) 

[c p 2 - (w 2 /k 2 )] [ Cj 2 - (w 2 /k 2 )] [c 2 l - (w 2 /k 2 )] = 


? it 2 2 

w r 2 , 2 .. 2. , 2 C A 

[Cs - (u/k )] Cj —z- 

* 4 


(2.3) 
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where w is the wave frequency, k the wave number, and w/k the phase 
velocity. c p , Cj, and c gL are the fast, intermediate, and slow hydro- 
magnetic wave speeds, defined by 

- 1 1/2 


■ 


C A + C S 


F 

C SL 

2 2 2 d 
C I = C A C0S ® 


2 + 


2 2 
C S + C A 


2 2 2 
- c„ c. cos 0 
S A 


(2.4) 


0 = cos * k j; /k is the angle between the wave vector and the magnetic 
field, c^ is the Alfven speed 

C A 2 = B o 2 /4ttNM + (2.5) 

where N. is the equilibrium number density, c g , the sound speed 


. 2 n T+ + V.T' 

C S " M 


( 2 . 6 ) 


4 » 

and y+ s the ion or electron ratio of specific heats. 

2 2 2 

The term k c^ represents finite ion inertia. The current 

velocities in hydroma gnetic waves are the order of c^ ; for wavelengths 
comparable with the ion cyclotron radius based upon c., ion currents decouple 

A 

from the waves. Since ion inertia is decoupled from moving lines of 
force, the effective Alfven wave speed increases; since the degree of 
decoupling is wavelength dependent, the propagation is dispersive. When 
k increases to w /c , finite electron inertia reduces the electron 

P_ 

current, thereby slowing the effective Alfven speed. Since this occurs 
at wavelengths much shorter than those of interest here, finite electron 
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inertia terms have been dropped in (2.2). When kc^/Q^ 0 , (2.2) 
reproduces the single fluid hydromagnetic result of three uncoupled 
nondispersive waves; finite ion inertia couples and makes dispersive the 
three hydromagnetic waves. 

Before considering finite ion inertia specifically, we compare the 

drift velocity with single fluid hydromagnetic wave speeds. For example, 

2 2 

when Cg / is small, the fast wave speed is given approximately by 

2 2 2 2 2 ± 
to = k + Ki. c g , where kj_ = k sin0 . Since w/k^ » v^ (when 

« 1 ) for all angles of propagation, the fast wave weakly couples 

to drift motions. In general, hydromagnetic velocities are comparable 

with the thermal speeds, and so drift effects may be neglected. Propagation 

almost perpendicular to the magnetic field (k^» k (i ) is an exception. 

For example, io = k tl c A for the intermediate wave for all angles of 

+ 

propagation; by choosing k i( /k^ - R + /L « 1 , we can reduce w/k^ 


enough to be comparable with v 


A similar argument applies to the 


2 2 2 2 2 

slow wave; for example, when c g /c^ « 1 , to = k t; c g ; again to is 

independent of kj.. , so that increasing k^ decreases to/k_ . Thus, 

the intermediate and slow waves have drift modifications when their 

perpendicular wavelength is much less than their parallel wavelength. 

The above argument suggests that we need consider the effects of 

finite ion inertia only upon the nearly perpendicular intermediate and 

2 2 2 

slow waves. When k^ » k (< , Cp » to /k for these waves, and the 
bi-cubic (2,3) reduces to a biquadratic, removing the fast wave from 
consideration. 


[Cj 2 - (w 2 /k 2 )] [c 2 l - (w 2 /k 2 )] = [c s 2 - (u 2 /k 2 )] r 

1c 


( 2 . 8 ) 



- 8 - 


.22 2 
k C A C I 

where r = — -- ^ is a finite ion inertia parameter; r / 0 couples 

Q c„ 

+ F 

the intermediate and slow branches. The sign of the dispersive coupling 
depends upon the ratio of to/k and Cg . Solving (2.8), 


w 


2 2 2 2 2 2 

1 ) C I + C SL * rC S A C I “ C SL + rC S 

+ x 


1+r 


When 0 < r « 1 , the intermediate wave is given by 


.2,2 2 

"s 

1/2 

. . 4c SL tc S - C I )r 


1 + 2 2 2 
[>=1 »♦*) - 

. 


(2.9) 


(0 


~2 " C I 


1 + r 


'I 


c 2 _ c 2 
I C SL 


( 2 . 10 ) 


Since c^ - Cg^ > 0 always, intermediate wave speed increases as r 


. 2 2 2 2 
increases if Cg > c^ , and decreases if Cg < c^ 


2 _ 2 2 fl 
c T = c. cos 0 
I A 


2 2 2 2 2 
implies that when Cg /c^ < 1 , waves with 0 £ cos 0 < c g /c^ speed 

2 2 2 

up, and those with Cg /c^ £ cos 0 slow down with increasing r . 
Similarly, when r is small, the slow wave is 


(0 


2 = C SL 


1 - r 


2 2 

C S “ C SL 
2 2 

C I - C SL 


( 2 . 11 ) 


2 2 

Since Cg > c gL , the slow wave decreases its phase velocity with 
increasing r . As r •+ 00 , the intermediate wave approaches an 
isotropic sound wave 

2 2 2 
w /k = c„ 


( 2 . 12 ) 
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and the slow wave approaches a wave resonance (w cos0) 


since 



2 2 
C F C SL 


2 2 
Cp c SL 

2 2 2 , 2 

(k c A /fl + ) c g 

2 2 2 q 
c. C c cos 0 . 
A S 



cos 2 0 



(2 . 13) 


2.3) Polarizations 

Here, we estimate the parallel wave electric field 6E t; , and the 
compressional component of the wave magnetic field, 6B >o . The ratio of 
6E u to the electric field component in the direction of the perpendicular 
wave vector, 6E X , is (Stringer, 1963; Formisano and Kennel, 1969) 



T [1 - (k 2 Cj 2 /w 2 )] cos0 sin0 

[(w 2 /k 2 Cg 2 ) - 1] + t sin 2 0 [1 - (k 2 Cl 2 /u 2 )] 


(2.14) 


y_t 

Y + T + + Yj“ 

2 2 

where the appropriate u) /k from (2.3) must be substituted in (2.14). 

In the hydromagnetic limit (kc^/fi^ = 0) , u = kCj precisely, 
whereupon (2.14) implies <5 e,j = 0 ; this means hydromagnetic waves have 
no Landau resonance interactions. Thus, while their propagation may be 
affected by drift terms, hydromagnetic Alfven waves can have no drift 
instabilities (Kennel and Greene, 1966). Substituting (2.10) into (2.14), 
6E./6E, is, to lowest order in kc & /£l , 

II -A. *3* 
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«E h 


.2 2 
k C A 


(2 



2,2 2 
C F (C I ' C SL^ 


T cosQ sin0 


(2.15) 


suggesting that strong resonant Alfven wave-particle couplings occur when 
kc^/S 1 + - 1 . By contrast, the slow wave has SE m j* 0 , and consequently 
Landau resonance interactions, even when kc^/ft + « 1 , For example, 
when Cg /c^ « 1 , 6E t /6E^ - tan0 = k/k x , implying that the slow ion 
acoustic wave is electrostatic, and has a zero wave magnetic field. 


In this section, we removed the magnetosonic wave from consideration 
2 2 2 

by assuming c p » w /k ; in the kinetic theory calculation, we will 
do so by constraining the wave compressional magnetic field amplitude. 


6B , 


to be zero. (In hydromagnetics, the fast wave has a nonzero 5B 


i, » 


2 2 

while the intermediate wave always has <SB,j = 0 , and when c g /c « 1 
the slow wave has a small 6B„. . Thus, assuming 6B ti = 0 picks out 
only slow and intermediate waves.) If k_^ is in the y-direction, the 
hydro magnetic intermediate wave has SB^ f 0 . From two-fluid theory, 


6B “ , k2c i 2 % 

= (i — 1_ j 

X b) 


0) 


,2 2 

k C A 


fi + tan0 
u> 


(2.16) 


substituting (2.10) into (2.16), assuming cos0 « 1 . 


6B .; 

6B 

x 



(2.17) 


Thus, both (2.16) and (2.17) indicate the hydromagnetic Alfven wave has 

2 2 

6B, = 0 . When c c /c. « 1 , <5B is small even when kc./ft - 1 • 

ti S A ' it A + 

2 2 

The slow wave is electrostatic when Cg /c^ « 1 , implying a small 6B u . 

. We may argue alternatively as follows: For perpendicular propagation, the 

slow wave reduces to a discontinuity across which pressure is conserved. If 
2 2 

C S /c A is small, the fractional change in particle pressure across the slow 
discontinuity is larger than that in the field strength. 
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2.4) Summary 


1.) Finite ion inertia dispersion causes the intermediate wave speed to 

‘i 

approach c g at short wavelengths. If c g > Cj , the phase speed 

increases to c~ with increasing kc./ft . When c T > c„ . the 
o A + IS 

phase speed decreases to c c with increasing kc./Q . The first 
case is appropriate to the magnetosphere, where particle pressures 
are large. Finite ion inertia decreases the slow speed; the slow 
wave frequency is alv/ays less than cos8 . These features are 
illustrated in figure 1. 


2.) When k /k = R /L « 1 , intermediate and slow waves can interact 

*i + p 

with drift motions. Magnetosonic waves have no significant drift 
interaction at any propagation angle. 


3.) The slow wave has a nonzero 5E., , which permits Landau resonance 
interactions . 


4. ) Drift Alfven instabilities could explain magnetically polarized 

micropulsations . However, <5E 1( = 0 in the hydromagnetic limit, 
implying hydromagnetic Alfven waves have no Landau coupling to particles. 

5. ) Finite ion inertia couples the slow and intermediate branches, causing 

an exchange of polarizations. With increasing , the inter- 

mediate Alfven wave acquires a dE^ (from the almost electrostatic 
slow wave), permitting Landau interactions at short wavelengths, 

kc A /!! . * 1 • 
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6.) 5B f , is large for the fast wave, zero for the hydromagnet ic 

intermediate wave, and small for the slow ion acoustic wave when 
2 2 

c„ /c. « 1 . These relations hold for finite kc JO. so lone 

S A A + 

2 2 

as Cg /c^ is not too large. 
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3 . Kinetic Theory of Drift Wa v es in a Uniform Magnetic F iel d 

3.1) Introduction 

Two-fluid theory, even extended to spatially inhomogeneous plasmas, 
describes neither resonant wave-particle interactions, nor perpendicular 
wavelengths near the ion thermal cyclotron radius. Therefore, we turn to 
kinetic theory, making as many simplifying assumptions as possible before 
performing the calculation; these are described in section 3.2. In 
section 3.3, we derive a dispersion relation for coupled drift ion 
acoustic (slow) and drift Alfven (intermediate) leaves. Section 3.4 
closes with a comparison with two-fluid theory and with previous work on 
drift instabilities. 


3.2) Approximations 

a.) Straightline magnetic geometry. 

We take the equilibrium magnetic field in the z-direction of 
an (x,y,z) Cartesian coordinate system. The equilibrium is uniform in 
y with gradients only along x . (For radial magnet osph eric gradients, 
x 'v radius, y 'v longitude, z direction along magnetic field.) By 
neglecting field curvature, we cannot find the parallel mode structure, 
important for relating ground and satellite micropulsation experiments. 
[Cummings et al. (1969) have in fact calculated the mode structure of 
hydromagnetic standing Alfven waves in the period range of interest here.] 
However, the instability is probably qualitatively described by straightline 
geometry. In a paper on electrostatic drift waves in curved magnetic 
geometry, Rutherford and' Frieman (1968) indicate that curvature introduces 
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no new electrostatic inodes or instabilities but only modifies details of 

those previously found in straightline' approximation. No curved field 

calculations to date (e.g., Rosenbluth, 1968; Rutherford and Frieman, 1968; 

Liu, 1969) treats either the electromagnetic polarizations or perpendicular 

wavelengths near the ion cyclotron radius, important here. Nonetheless, 

from general geometry experience, we can extract two points. First, 

Landau interactions in straightline geometry correspond to bounce resonance 

interactions in general geometry. Since the most unstable drift waves have 

long wavelengths parallel to the magnetic field, whose upper limit is the 

length of the line of force, the k (< v (l in the Landau resonance condition, 

u - k.v (1 = 0 , corresponds approximately to the bounce frequency of a 

particle mirroring near the earth. Second, general geometry calculations 

involve quantities averaged over the length of the tube of force, i.e., 
f dll 

<Q> = j -j-g-j — Q , where dll is an element of length along the line of 
force. This averaging process emphasizes those points where } B | is small, 
suggesting that equatorial plane parameters make the appropriate correspondence 
between slab geometry calculations and the magnetosphere. 

b.) Neglect of ionospheric attenuation and coupling. 

We neglect wave energy losses due to resistive ionospheric 
attenuation or to coupling to the ionosopheric wave guide. To overcome any 
such losses and ensure instability, the equatorial plane growth rate must 
be finite and positive. However, if ionospheric losses are not strongly 
frequency dependent, the equatorial growth rate determines which modes 


are most unstable. 
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c.) Assumption of small magnetic inhomogeneity. 


We assume 


+ i d 0 -ti ■n H 

v n « , This requires that — d3c~ >:> “ 


Jin B 


dx 


The hydromagnetic pressure equilibrium condition, V[p + (B /87O ] = 0 


implies for slab geometry, that 


d Jin 
dx 


d Jin B 
0 dx 


, where 


3 = 8iTp/B . Thus, if 0 « 1 , |v D + | « |y *1 . For the inner edge of 

3 

the electron plasma sheet, where N ^ 1/cm , T = 3 keV , B = 60y , 

0-0.3 , suggesting that jv^| =0(5) | v D | . In other words, when the 

electron 0 < 1 , the logarithmic gradient of the magnetic field is roughly 

that of the dipole field. At 1 = 6, the dipole scale length is = 2 , 

1 

whereas the electron scalelength is ~ R £ (Vasyliunas, 1968a). Our 
analysis is restricted to such sharp particle distribution gradients. 


d. ) Annihilation of fast magnetosonic wave. 

The magnetosonic mode is always faster than typical drift speeds, 
and does not have large drift modifications; its presence only clutters 
drift calculations with irrelevant terms. We remove it in a dvance, by 
constraining the modes we do treat to have at most rotational magnetic 
field components, i.e., by requiring 6B (j = 0 . Recall that only the fast 
wave has a significant nonzero SB t| . An alternative argument is that the 
fast v/ave adjusts the total pressure time scales shorter than drift time 
scales, leaving only small 6B ([ adjustments to the other two waves. 

e. ) Localization of spatial eigenfunctions. 

Because the plasma equilibrium is inhomogeneous, an 

integrodifferential equation describes, in general, the x-dependence of 
the waves' spatial eigenfunctions. For short x-wavelengths relative to 
equilibrium scale lengths, the WKB approximation, wherein the wave fields 
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vary as exp 1 


i / X k (x 1 ) dx' + k y + k z 
x J y J z 


artTJ 


may be used. Since the 


instability energy source is a localized spatial gradient, the fastest 
growing modes should have eigenfunctions localized near the maximum plasma 
gradient. For these waves, the WKB dispersion relation has the form 


x_ 

r 2 


dx' k x (x , ;w\k y ,k z ) = nir 


W = W " 0 


(3.1) 


where the wave amplitude must tend to zero as x ± °° . Rukhadze and Silin 
(1964), Mikhailovsky (1967), Krall (1968), and others have argued that a 
good estimate of the eigenfrequencies and growth rate can be obtained by 
examining the turning points, where k = 0 . Therefore, we set k = 0 

X X 

in advance, creating an algebraic dispersion relation, rather than an 

integrodifferential equation with boundary conditions. 

Combining the SB,. = 0 and k =0 conditions with Faradav's law 
° a x 

implies that SF.^ , which carries the fast wave, is zero. The only 

remaining wave components are 6E , 6E E 5E , SB , k = k x , 

y z x y 

and k = k . 

Z M 


3 . 3) Calculation of Disp e rsion Relation 

+ 

The equilibrium distribution functions F q ~ depend only upon the 
particle constants of motion. Time invariance and translational 
invariance along B imply the particle energy and parallel momentum, 
respectively, are conservedjin nonrelativistic plasmas, v^ and v„ 
are equivalent constants of motion. From y-invariance, p^, - M + v^ + eA^/c 
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is conserved. For a constant magnetic field, the vector potential 

Ay = 1 B ! x , so that X = x + (Vy/fy.) is an equivalent constant of motion. 

We choose F q to be locally bi-Maxwe Ilian: 



r m. 


r m, 1 

,v ) - N 

± 


± 

4 * 


+ 


. 2 itT x”. 


L 2Tr vJ 


ll/2 


exp - 


M +V 


2T ; 


M + v 


2T, 


2 T 


(3.2) 


+ — 

N(X) is the electron and ion number density, and T “(X) . and T. (X) , 

the electron or ion perpendicular and parallel temperature, respectively. 
The equilibrium particle orbits in the frame for which equilibrium electric 
fields are zero, are: 


v x = v j. cosOM+t) J v y = Vjl sin(i|»-n ± t) , v z = v 4j (3.3) 

v x 

x(t) = x Q - [sin(<J>-fl ± t) - sim{>] 

+ 

v j. 

y(t) = y o + if [cos Ol»-n ± t) - sinfj (3.4) 

z(t) = z o + V w t 

where if) is the Larmor phase at t = 0 , and (x o ,y Q ,z o ) is the 
particle location at t = 0 . Finally, 

Vj sin(ty-fi + t) 

X = x(t) + — *5 —— (3.5) 

4 x 

i fk®x-wt) 

For small amplitude harmonic wave fields, varying as e - - , 

the wave electric field satisfies 

2 

k x (k x 6 E) + ~ e • 6 E * 0 (3.6) 

•w ++ Z •« *!** 

C - 

where the dielectric tensor e is defined by (Stix, 1962) 
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d\ v(6f + - 6f~) (3.7) 

Sj is the perturbation current , and Sf , the perturbed distribution 

. \ 

function. We treat one species of singly charged ions. The dispersion 
relation involves finding 6f~ , multiplying by v , integrating over 
velocity space, expressing 6J in terms of 6E , and solving the 
determinant of (3.6). 

6f“ can be found by the method of characteristics (Rosenbluth, 

1965; Mikhailovsky, 1967), wherein 6f~ at t = 0 is given by an 
integration following unperturbed particle trajectories: 


e*5E » 6E + 


4-rri 

to 


6j * 6E + 


4tt1 

CO 


I 


6f [x(t=0) ,v(t=0) ,0] 


0 



-CO 


v(t') x SB 
6e * - — — 

c 



exp (i tk_^y (t ' ) + k >j z(t*) - tot']) (3.8) 


where t' is a time index following particle trajectories. Im <o = y 

is assumed positive to ensure convergence of the integral at t" = - °°. 

Assuming SE^ = SB^ = SB^ = 0 , using Faraday's law to express <$B^ 

in terms of SE and SE , and making the change of coordinates 
y z 

+ + ^ 

F q [x(t),v(t)] -*■ F q (X,v x ,v m ) so that F q is a constant in the t' integration: 


M. 


v x SB 


6E + 


3F 


0 = M s.in(^Sl + t) + N 


3v 


where 


= + 


3F “ (k.,E - k.E ) 

O V « V x z 






y 

CO 


3F 


v. 


H W 


3v 


it - e 

N = + ft 


+ 4 * 

3F E 3F " 

° y _ o 


n 3v 


n 
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For the locally bi -Maxwellian distributions, (3,2) 


M~ = ± 


v, 

e -*■ 


M± Cc^ ± ) 2 
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M + ' 
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k v ' , 
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(3.11) 
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where the thermal speeds c^ + , c jj+ are 


<x± - CT 1 ± /M ± ) 1/2 ■ 


s± ■ CT^/M ,) 1 / 2 


(3.11a) 


x t* * 4 * " 4 * 

and the thermal anisotropy A~ is (T,, - - T 1 _ - )/T u “ . Note that 


3F 


o 3 En N 


3X 3X 


l - n, 


. 3 An T. 

X JL. 

* 4 . = 3 in N 


1 , 1 n. 

1 " 2 t v j/ c x±) ~ 1 
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(C J 


4,± 

(3.12) 


a 3 in N 


Using (3.4), we may write 


ik'xft'j-iot ’ xk*x 

„ ~ ~ o 

e = e 


T J J e 
L n m 


i(n-m) [i|h- (it/ 2)] -i(X+nft + )t' 
e 


(3.13) 


X s to - k ft v ( . 


d = J (kjV./n.) , m,n any integers 

li SI — 

where we used the Bessel function identity e* a S * n< ^ = £ J n (a) e * n< ^ 
twice. Combining (3.13) with (3.9), we perform the orbit integration 
noting that F q _ (X,v^,v) is independent of t' : 
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The perturbation currents are 
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(3.16) 


Equations (3.15) and (3.16) contain v^-integrations over F q » which 


contains directly and also indirectly through X . In evaluating 

(3.15) and (3.16) we assume the equilibrium scale length is much longer 


+ + 

than the ion cyclotron radius, so that F q (X) - F q (x) + 0(v^/fi + L j , 


and 3/3X - 3/3x , thereby neglecting the indirect v. dependence. 

When w , y , and the "bounce" frequencies k v tt of typical particles 
are well below the ion cyclotron frequency ^ X/£2 + « 1 , and we may expand 
the summation of cyclotron resonance denominators : 
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n X + nfi + X 
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(3.17) 
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where 6 =1 if n = 0 and zero otherwise. Equation (3.17) implies 

that we neglect all cyclotron resonance interactions. For perpendicular 

wavelengths near the ion cyclotron radius, none of the ion Bessel functions 
are small; thus we must sum contributions from all n . However, this is 
faciliated using 


l nJ 2 (1 ~— J 

% n ^n 


n 


n 


= 0 ; 


I J 2 = 1 ; 
*-n n 


I (l- 5 ) J : 

L n v no'' n 


1 - 



Keeping the first nonvanishing terms in A/fi + , (3.15) and (3.16) reduce to 


6J 


= i I (±e) f d 3 v V J 
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(3.19) 


The corrections to (3.18) and (3,19) are at least 0(A/ft + ) smaller than the 
terms kept. Note that the y-current, which is responsible for Alfven waves, 
has no 1/A Landau resonance terms, whereas the z-current, which leads 
to slow ion sound waves, does. 

The relevant components of e are: 

/ * 
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(3.20c) 
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G , which arises from v^- integral ion of (l-J o )F q ~ is 


+ 1 - e‘ z I (z") 1 - r 


G" = 


(3.21) 


where z~ = k ^T,~/M.R ^ , T~ - e Z I (z ) , and I is a modified 
-*• ■** ± ± n n n 

4* + 

Bessel function of the first kind. As z~ •* 0 , G -*■ 1 ; for 

+ + ± 
z~ » 1 , G" 'v 1/z 

w* + is a drift frequency: 


ur 
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wfi. 
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(3,22) 


x x „ X 

where the arguments z of i Q , have been suppressed. Dg(co,k) = 0 
is the electrostatic drift wave dispersion relation; 
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D(,(w»k) = £ - 9 — 


+ >" c i; ± 
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2 T 2 


d v v. J Q * [F o ~ - (k^c; ± /^ ± )(3F o -/aX)] 


k.w N 

it 


(w/k.) - v, 


(3.23) 


Since, in reaching (3.20), we set all odd v ;i moments of F q ~ zero, 
we cannot treat instabilities from currents and heat flowing along lines 
of force. 

1/2 

Dropping terms small in M_/M + (the largest of which is 0(M_/M + ) 

comparable with terms of 0 (A/Q^) already neglected), using 
2 2 2 2 

fi + /u> = c^ /c and dropping the unity terms in (3.20a) and (3.20c), 

2 2 

which are of order c A /c « 1 , we reduce (3.20) further: 
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r • V“- k > 

4 * 


(3.24a) 


(3.24b) 


(3.24c) 


n i> i + 

AT '£ 2 G (Tj - TJ") is an effective temperature anisotropy. 

To correspond to frequencies between the ion and electron bounce 
frequencies, we assume k c |t+ « 0 ) « k^c^ . For the intermediate speed 
to be smaller than the electron thermal speed (c A < c {) _) , 8 ttNT ( ,"/B 2 = 

3 > 2(M /M + ) . (Thus, overall, we require 1 » 3 ~ » 2(M /M + ) .) 

For the slow speed to be larger than the ion thermal speed (Cg > c + ) , 

„ 4 . 

T |( /T )( > 1 . If, moreover, we assume y/u « 1 , then, for ions 


w n ~ — ° lr » » -ir- « •••’ - ( 3 - 25a > 


and, for electrons 


— -p . — ^ = - — (1 + — + . . .) - iiT signk 6[(w/k ) - v ] 

(co/k,,) - v, v (j k w v (( s " " 

where the Dirac delta functions signify Landau resonance wave-particle 
interactions. Using (3.25), Dg(w,k) becomes 
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w k n C S 
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Unlike the fluid definition (2.9), here c c = T /1M . u is an ion 

S *» + + 


(3.25 


(3,26 


drift frequency: 
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w is the corresponding electron drift frequency 
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In section 4, we will use the fact that when z" « 1 , to is a pure 
number density drift. 6 contains resonant ion and electron effects: 
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Using (3.20) - (3.29), we may find the dispersion relation: 


[G + u 2 - uxi>* + - k.,^ £ *] [r o 'a)^ - S_w - k (J z c s^( r D ^ ” + ^ 


2-2, rr -.2 


2 _ 2f„ * 


----- - k.. 2 c 2 [gV - ww* ] (3.30) 

a 2 " s 


where c^ is an effective Alfven speed reflecting modifications by thermal 
anisotropy ("firehose") and finite cyclotron radius effects, 

d A 2 ’ C A 2 " " /M . • 
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3.4) Discussion 

Let us compare (3.30) with homogeneous two-fluid theory results. 
Neglecting thermal anisotropies and spatial derivatives, (3.30) reduces to 


,2 2 

r 2 ,, 2 2 rr- 2 .2 2 „ + K x C A , 2 2 2 

t“ - ( k „ C A /G )] [r o“ - k u C S r o * lS1 - — k u c s “ 


(3.31) 


f. t + + + 

Neglecting o , and noting that when z~ « 1 , G and F approach 

o 

unity, (3.32) may be rewritten 


( 2 n 2v n (0 

(U /k ) ] * “2 c s 
k 


. 2.2 2 
k k , c. 

2 « A 


using c| L 2 c g 2 (k it 2 /k 2 ) when Cg 2 /c ^ <<: 1 • When » c g 2 » u 2 /k 2 , 

(3.32) and (2.8) reduce to approximately the same expression, with 

kj_ 2 c^ 2 /n + 2 in (3.32) replacing k 2 c^ 2 /0 + in (2.8). Since k^/k 2 1 , this 

difference is unimportant. Thus, kinetic theory reproduces homogeneous 

2 2 + 

two fluid theory results when Cg /c^ and z are small, and resonant 
particles may be neglected. 

Finite cyclotron radius (FCR) effects, not described in fluid theory, 
may be evaluated schematically by neglecting finite ion inertia in (3-31), 
whereupon when z « 1 , the intermediate and slow waves are given by 
w 2 a k t 2 c^ 2 /G + and w 2 * k 54 2 Cg 2 r Q + respectively. Both FCR and finite ion 
inertia increase the intermediate speed and decrease the slow speed as 
k. increases, since both G and T decrease as z increases. 

6 represents Landau wave-particle interactions, which in homogeneous 
Maxwellian plasmas can only damp waves. Spatial inhomogeneity, however. 
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permits instability in certain cases. It also modifies the real part of 

the wave frequency, through to* and o>, . 

* » — 

(3.30) reduces to the results of Mikhailovsky and Rudakov 

4 

(1963) , who kept finite inertia coupling but assumed z « 1 and to 

those of Kennel and Greene (1966) who treated hydroraagnetic drift 

+ 2 2 
instabilities (z - kc A /& + - 0) for arbitrary c g /c^ . When z = 

ri _ = ri"' = AT = 0 , (3.30) reproduces the results of Mikhailovsky (1967). 
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4 • Drift Instabilities of the Auroral Electron Bou nd ary 


4.1) Alfven Thermal Gradient Instability 

The auroral boundary electron temperature gradient is much larger than 
the number density gradient (Vasyliunas, 1968a). The ions have no sharp 
gradient near the electron boundary (Frank, 1968). Thus, we keep only 
the electron temperature gradient. Since when z « 1 , to - 0 for 
small number density gradient, (3,31) reduces to (3.32), the homogeneous 
plasma dispersion relation, except that <5 has thermal gradient contribu- 
tions. When y/h) « 1 , (3.32) reduces to separate equations for y and w : 



f“ 2 - k 2 c s 2 Vi 


2 , 2 2 ,, 2 2 /r . 2 . 
“ C S C A /fi + 3 


(4,1a) 



2 ,,2 2 ,_+ , 
w - (kj, c A /G ) 

r 4 ,, 4 5 2 r + . +. •, 

[ “ - (k i C A C S P o /G 51 


(4.1b) 


where AT is assumed zero, y > 0 signifies instability. In order that 
2 2 2 + 

y ? 0 , to i lc^c^ /G ; this condition is always satisfied for the slow 
wave, and for the intermediate wave only when finite ion inertia is 
accounted for. 

+ 2 2 2 
When z > 0 , both FCR and finite ion inertia ensure w > k c^ 

2 2 2 

and to < k u Cg for the intermediate and slow waves respectively. Since, 
2 2 

moreover, Cg / c A < 1 , we may decouple the two waves in (4.1a) by 

2 2 2 + 

assuming w » k Cg r Q for the Alfven wave, whereupon 
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2 2 2 

Similarly we may assume to /k u c + » 1 , permitting neglect of ion 
Landau interactions. Then, substituting (4.2) into (4.1b) 


Y 

tO 


- /tt/ 8 


0) 


exp[-w^/2k ll ^c 

•T,. 



A r i , 2 n 2 2. , 

— [1 - (to /k h c_ )] 


(4.3) 


A is the electron temperature drift frequency 


k jL c - d An T - _ T" /T C S 

2 dx “ T + 2L T 


(4.4) 


where L^ = (d An T /dx) * . IVhen A = 0 , the Alfven wave is damped. 
A sufficient condition for instability is 


1 + ~ [1 - (w 2 /k ,fc_ 2 )] < 0 (4.5) 

2 2 2 j 

Since by assumption to /k,j c < 1 , A/to < -1 is a necessary condition 

for instability, implying that the perpendicular phase velocity to/k^ and 

2 d £n T” 

the thermal gradient drift velocity c /2Q, — are oppositely 

directed. } A/to | > 1 implies 



Large leads to instability. T /T + large, c g /c A lar 8 e » an< ^ a lar 8 e 

ratio of parallel wavelength to thermal gradient scalelength (k ( Lj, << 1) 
all contribute to instability. When z + « 1 , - 1 - z* and 

(1 - T ) - l/z , so that instability does not occur for z 0 . 

z + must surpass a threshold value for instability; for instance, when 



-29 


A. = 2 , and T~/T + - 2 , z* > 0.8 for instability. Thus, except for 

A 

very large , instability sets in at wavelengths near the ion cyclotron 


radius . Since when z » 1 , both A and to grow as 


Js 


(for k 


constant) , a necessary condition for large z and thus any z to be 


unstable is obtained by setting T* -► 0 in (4.6), A. > f\ + T /T + 

O A 

4 » 

For k t - constant, y increases with increasing <o and z until 
to * k (j c_ , when y changes sign. When k t| corresponds to the lowest 
parallel mode, the unstable wave frequencies are below the electron bounce 
frequency and above the larger of the Alfven bounce frequency, k^ c^ , or 
the proton bounce frequency, kjjC + (since ion Landau damping would kill 
the instability). Since y maximizes just before it changes sign, the 
most unstable frequency is near but below the electron bounce frequency. 

In curved magnetic geometry, the growth rate peaks even more strongly 
near but below the bounce frequency since w/]k )t | c is replaced by a 

3 

term of 0(u)/k #i c ) (Rutherford and Frieman, 1968). 


4 . 2) Ion Sound Thermal Gradient Instability 
2 2 2 

Assuming to « k (| c^ , the slow ion sound wave frequency and 
growth rate are given by 


2 

to 


.2 2 _ + 

k . c s r o 

2 2 2 *T 

1 + (k/cgX ) G 





i + (T~/T + ) (1 - r D + ) 


(4.7) 
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(T"/T + ) 3/2 
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(4.8) 



- 30 - 


••• 4 * 

where we assumed w/k^c very small, true except when T /T is 
unrealistically large. A necessary condition for instability is A/w < -1 , 
which requires 



" 1 1/2 


2k « L T z + [l + (T"/T + )(l - r + )]_J 


( 4 . 9 ) 


This instability criterion, while less stringent than that for Alfven 
waves, still requires that z* be finite. Therefore, the unstable 
frequencies lie below k |( Cg . As z increases, w/k^ decreases to c^ , 
where the basic assumptions of this calculation break down. Thus, the 
lowest parallel mode will probably be unstable only above the ion bounce 
frequency . 
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5. Summary 

Only the intermediate and slow wave branches have strong drift 

interactions. To do so, they must propagate nearly perpendicular to the 

2 2 

magnetic field (k^/k^ - R + /L^ « 1). We restrict ourselves to c g /c^ 
small. Then, the slow wave is the ion sound wave, which always has a 
parallel electric field, and therefore, Landau interactions. On the 
Alfven branch, 6E., f 0 , permitting Landau resonance instabilities^ only 
at short perpendicular wavelengths, where finite ion inertia coupling 
permits an exchange of polarization with the ion acoustic wave. 

We derived the properties of drift Alfven and ion sound waves, 
coupled by finite ion inertia, in section 3. The various limitations of 
this theory are outlined in section 3.2; in particular, we neglect 
magnetic curvature, and the effects of magnetic gradient drifts. The 
second assumption restricts us to very sharp distribution function 
gradients, such as the auroral electron boundary. In the remainder of 
the magnetosphere, the density and magnetic drifts are comparable. In 
this case, the appropriate drift calculations must be generalizations of 
Rosenbluth, Krall and Rostoker (1962), who treated k ( | = 0 only, or of 
Kennel and Greene (1966), who treated only z*<< 1 . Mikhailovsky and 
Fridman (1968), have recently calculated the properties of high 3 drift 
waves with V D ^ 0 . 

A great simplification occurs for the auroral electron boundary, 
since we may neglect all spatial gradients but that of the electron 
temperature. Here, the drift terms disappear from the real part of the 
dispersion relation, leaving it very similar to homogeneous two-fluid 
theory; however, the growth rate has a destabilizing contribution from 
the temperature gradient. Let us now compare the qualitative predictions 
of this theory with experiment. 
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5.1) Satisfaction of Instability Criterion 

We compute A A for parameters appropriate to the postmidnight auroral 

zone equatorial plane. For T~ = 4 keV , T + = 1 keV , B = 70 y » 

3 7 a 

N - 1/cm , Cg = 8.6 x 10 cm/sec , - 4.8 * 10 cm/sec , and 

Cg/c^ = .18 . The longest parallel wavelength available is the length of 

the line of force, “ 20 at L = 6 . For this lowest parallel mode 

A a - 1.1 Rg/L T . From (4.6), A A > s/l+ T /T + = 2.2 implies instability. 

Thus, if L t < 0.5 Rj, , Alfven wave instability is possible. Observed 

values of L^, fall within this range (Vasyliunas, 1968a). On the other 

hand, higher parallel modes are not likely to be unstable for this range 

of parameters. 

5. 2) Unstable Frequencies and Wavelengths 

When only the lowest parallel mode is excited, the unstable frequency 
range should be narrow, lying above the Alfven bounce frequency (k|, c A ) 
and the thermal electron bounce frequency (k ( ,c ) . For the parameters 
chosen above, these two outer limits correspond to periods of 4 and 25 
seconds. Since the growth rate maximizes near but below k (V c , waves with 
- 10 second periods should dominate the spectrum. McPherron et al. (1968) 
have characterized the observed micropulsations as "band-limited" in the 
5-15 second period range. 

Corresponding to the narrow unstable frequency range, a small range 
of perpendicular wavelengths, approximately 1 < 77 < c /c A (cf. Section 

4.2), should be unstable: perpendicular wavelengths near the thermal 

proton cyclotron radius in the equatorial plane, some tens of kilometers. 
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5.3) Location, Localization, and Propagation 

The association with the inner edge of the plasma sheet suggests that 
Alfven micropulsations will be located equatorward of the main auroral 
activity, in agreement with the location of auroral light pulsations 
(Cresswell and Davis, 1966). The fastest growing modes will be localized in 
radius to the region of maximum electron temperature gradient. If the 
equatorial thermal scalelength is -j \ ~ km , the ionospheric 

scalelength is the order of 100 km, consistent with the localization of 
energetic electron precipitation pulsations (Barcus et al., 1965; Barcus 
et al., 1966; Parks et al., 1968a) and of 10 second auroral light modulations 
(Cresswell and Davis, 1966). 

Since the instability criterion (4.6) involves T'/T* , the steady 
convection electric field across the magnetosphere (Axford and Hines, 1961) 
could account for the observed diurnal occurrence pattern for 5-15 second 
micropulsations. Electrons are accelerated on the morning side and ions 
on the evening side (Brice, 1968). Thus, if 1/1 =1 in the magnetospheric 

tail, T /T >1 on the morning side and 1/1 <1 on the evening side. 

OB> *|* 

Choosing T - 1 keV and T - 4 keV with all other parameters the same as 
above, - .03 Rg/L^. , strongly suggesting that instability is not possible. 
5-15 second auroral micropulsations are observed primarily in the 0200-1000 
L.T. postmidnight sector (McPherron et al., 1968). Thus, while the auroral 
electron boundary is more pronounced on the evening side (Vasyliunas, 1968a) 
than on the morning side (Vasyliunas, 1968b), it appears that the electron 
to ion temperature ratio may play the controlling role in determining 
instability. 

If we assume the eigenoscillation is a standing wave in x (radius) 
and z (along the lines of force) , then energy propagates along y , or in 
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longitude. From (4.5), the unstable waves propagate opposite to the 
electron thermal gradient drift, or towards morning, with phase velocities 
less than A/kj_ - 10 km/sec, in the equatorial plane. This velocity, 
mapped onto the ionosphere, is - 1 km/sec. If the micropulsations modulate 
electron precipitation, then the phase surfaces of associated auroral light 
irregularities should progress eastward at near but below this velocity. 
Cresswell and Davis (1966) found pulsating auroras move eastward with 
0.1-1 km/sec velocities. The fact that no north-south motions were 
observed is consistent with a radially confined mode. 

5.4) Polarization 

Our calculation keeps only 6 B x , corresponding to a radial perturbation 
magnetic field in the equatorial plane, which maps into a north-south 
perturbation field in the ionosphere if there is no rotation of the 
polarization along the line of force. Accurate calculations of the x-dependence 
of the eigenfunction are needed to determine the actual SB^. , or longitudinal 
component. In any case, the rotational components of the perturbation field 
should be larger than the compressional components, though because c<Vc^ 
not always small, 6B vi is not necessarily strictly zero, as was assumed 
here. 5-15 second band limited micropulsations are observed magnetically 
on the ground (McPherron et al., 1968), consistent with a magnetic polariza- 
tion in space. 
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5.5) Estimate of Saturation Amplitude 

The development of the drift instability should smooth the electron 
temperature gradient driving it, by diffusing particles across lines of 
force. Assuming this process is incoherent, we can write an approximate 
heat flow equation 


3T~ _ „ 3 2 T' 

3t ~ U “ 2 
3x 


(5.1) 


where the heat conduction coefficient should be given dimensionally by the 

magnitude of the particle velocities in the wave (c6E /B q ) and the wave 

2 -1 

frequency D - (c5E /B Q ) u) . In the strong diffusion picture of the 

auroral electron boundary (Petschek and Kennel, 1966; Kennel, 1969; 

Vasyliunas, to be published) the minimum lifetime for electron precipitation 

T m is the scale tine for boundary formation. Assuming that thermal 

transport also has this scale time permits a crude estimate of the amplitude. 

2 3 

Equating T M = L T /D = 10 sec (for L = 6 and keV electrons), and 
noting <5E = (c^/c) SB , 


T M = (V C A )2 “WV " 2 


(5.2) 


8 

For L T = 0.5 Rg , w = 2ir/10 radians/sec , c^ = 2 x 10 cm/sec , 

2 

6B/B q - .03 , or 6B = 3 y when B q = 10 y . This wave amplitude ought 
to be dominantly rotational. Ground based measurements indicate 6B = 1 y 
(McPherron, private communication). Let us now ask whether Alfven micro- 
pulsations can create an electron precipitation pulsation. To achieve an 
effective congressional magnetic component the Alfven wave must carry 
particles into an increasing part of the dipole field. In this case, the 
effective compressional amplitude b , defined in Coroniti and Kennel (1969), 
is 
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b * 


3A 

R — 
O LR C 
£ 


(5.3) 


2 

where A^ is the perpendicular wavelength. Choosing B q = 10 y , 

7 

A - 30 km (approximately a 1 keV ion cyclotron radius), = 6 x io km, 
and L = 6 , we find b - 1.5 y , consistent with their estimate for 
electron precipitation pulsations with peak- to-val ley flux ratios of 2. 
Thus, if this picture is correct, whistler amplitudes ought occasionally 
to be modulated by Alfven micropulsations. Brody et al. (1969) have 
reported preliminary observations of whistler chords bursts, which repeat 
quasi-periodically with typically 6-11 second periods, in the postmidnight 
auroral equatorial plane. A schematic summary of these ideas is shown in 
figure 2. 


5 . 6) Ion Sound Thermal Gradient Instability 
According to (4.9), k^L^ need not be as small for ion sound 
instability as for Alfven instability, suggesting that if Alfven waves are 
excited, so will lower frequency ion sound waves. However, if Alfven waves 
are somewhat above marginal stability, their growth rate is larger than 
the ion sound growth rate. Then, if waves with the largest growth rate 
dominate the developed spectrum, or if ionospheric damping is important, 
the ion sound component may not be apparent. Since higher parallel ion 
sound modes could also be unstable, estimates of the unstable frequency 
spectrum must wait until curved magnetic field calculations with ionospheric 
damping are carried out. However an interesting, if extremely difficult, 
experimental question is raised: are there nearly electrostatic micro- 

pulsations? 

Before the question of whether drift instabilities can account for 
various different micropulsations can be answered, several theoretical and 
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experimental developments are necessary. Some of these necessary 
developments have been indicated by the present calculation. On the 
theoretical side, curved magnetic geometry calculations, including full 
electromagnetic polarizations, finite B, resonant particles, and short 
nonhydromagnetic wavelengths, must be carried out. The magnetic gradient 
and curvature guiding center drifts must also be retained, since the 
particle scalelengths will be longer than in the auroral electron 
boundary elsewhere in the magnetosphere. Even in the present calculation, 
neglect of the gradient drift presents some uncertainty, since the 
unstable waves propagate eastward with a velocity comparable with the 
guiding center drift velocity of Van Allen electrons of a few tens of keV 
energy. The loss of wave energy to the ionosphere should be evaluated. 

On the experimental side, the relationship of spatial gradients in the 
particle distributions to micropulsation activity needs to be elucidated. 

The question of electrostatic micropulsations is unresolved experimentally. 
Here the association with precipitation pulsations could prove quite 
useful, since particle modulations in the absence of magnetic oscillations 
could be evidence for electrostatic micropulsations. Nevertheless, 
despite these technical obstacles, drift instabilities appear to be a serious 
candidate for various micropulsations. An instability analysis of the 
auroral electron boundary, the one case where the particle gradient is 
experimentally well-defined and which fortuitously is open to the limited 
analytic techniques used here, has led us to propose a wave generation 
mechanism which is at least not inconsistent with the observations of period, 
spatial location and localization, and propagation speed and direction, of 
a common auroral micropulsation. 
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Fi gure Captions 

Figure 1. Oblique Hydromagnetic Waves' with Ion Inertia. 

These plots schematically summarize the effects of ion inertia 
on wave propagation nearly across the magnetic field in a homogeneous 
plasma. Figure la describes the Cg/Cj > * case ordinarily appropriate 
to the magnetosphere, and lb, the Cg/c^ < * case. In both cases, 
the fast wave has much larger phase and group velocities than the 
slow and intermediate wave. In both cases, finite ion inertia converts 
the intermediate wave into an isotropic sound wave (w/k = c g ) , and 
drives the slow wave to an oblique ion cyclotron resonance, w = fi cos8 . 
The slow and intermediate waves have a small compressional magnetic 
field component, 6B jt . EM denotes mixed electric and magnetic, and 
ES, approximately electrostatic polarization. 

Figure 2. Disturbed Auroral Electron Boundary. 

When the electron boundary is sharp enough, a drift Alfven 
instability creates - 10 second period micropulsations, primarily 
after midnight. When unstable whistler turbulence is already present, 
the micropulsations modulate the whistler instability, leading to 
finite amplitude precipitation pulsations observable in the aui’oral 
ionosphere , 



UNS 

ENT 






E HYDROMAGNE 
ION INERT 








-40- 


References 

Barcus, J. R. , and T. J. Rosenberg, Observations on the spatial structure 
of pulsating electron precipitation accompanying low frequency 
hydromagnetic disturbances in the auroral zone, J. Geophys. Res. , 70, 
1707, 1965. 

Barcus, J. R. , R. R. Brown, and T. J. Rosenberg, Spatial and temporal 

character of fast variations in auroral zone x-rays, J. Geophys. Res. , 
71., 125, 1966. 

Braginski, S. I., On the modes of plasma oscillations in a magnetic field. 
Spy. Phys. Dokl. , 2^, 345, 1957. 

Brice, N. M. , Magnetospheric and High Latitude Ionospheric Disturbance 
Phenomena. In: Physics of the Magnetosphere, R. L. Carovillano, 

J. F. McCloy, and H. R. Radoski, eds., D. Reidel, Dordrecht, 1968. 

Brody, K. I., R. E. Holzer, E. J. Smith, C. F. Kennel, R. W. Fredricks, and 
F. L. Scarf, Nightside Auroral Electric and Magnetic Noise in the 
Equatorial Plane, TRANS. Amer. Geophys. Union , 50 , 3, 1969. 

Coroniti, F. V., and C. F. Kennel, Electron precipitation pulsations, 

J , Geophys . Res . (this issue) . 

Cresswell, G. R. , and T. N ! . Davis, Observations on pulsating auroras, 

J. Geophys. Res, , 71 , 3155, 1966. 

Cummings, W. D. , R. J. Sullivan, and P. J. Coleman, Jr., Standing Alfven 
waves in the magnetosphere, J. Geophys. Res. , 74 , 778, 1969. 

D’Angelo, N., Role of the universal instability in auroral phenomena, 

J. Geophys. Res. , 74 , 909, 1969, 

Formisano, V. and C. F. Kennel, Small amplitude waves in high 3 plasmas, 

J. Plasma Phys. , 3_, 55, 1969. 



-41- 


Frank, L. A., On the distributions of low energy protons and electrons in 
the earth's magnetosphere. In: Earth's Particles and Fields, B. M. 

McCormac, ed.. Reinhold Book Corp., N. Y., 1968, p. 67. 

Kantrowitz, A. R., and H. E. Petschek, MHD characteristics and shock waves. 

In: Plasma Physics in Theory and Application, edited by W. B. Kunkel, 

McGraw Hill Book Co., N. Y., 1966, p. 148. 

Kennel, C. F. and J. M. Greene, Finite Larmor radius hydromagnetics, 

Ann, of Phys . , 38 , 63, 1966. 

Kennel, C. F. , Consequences of a magnetospheric plasma, Rev. Geophys. , l_ s 
379, 1969. 

Krall, N. A., Drift waves. In: Advances in Plasma Physics, edited by 

W. B. Thompson and A. Simon, Academic Press, N. Y., 1968. 

Liu, C. S., Interchange instabilities of the ring current belt, to be 
published, J . Geophys . Res . , 74, 1969. 

McPherron, R. L. , G. K. Parks, F. V. Coroniti, and S. H. Ward, Studies of 
the magnetospheric substorm, 2, Correlated magnetic micropulsations 
and electron precipitation occurring during auroral substorms, 

J . Geophys . Res . , 1697, 1968. 

Mikhailovskii, A. B., and L. I. Rudakov, The stability of a spatially 

inhomogeneous plasma in a magnetic field, Sov, Phys. JETP , 17 , 621, 1963. 

Mikhailovskii, A. B., Oscillations of an inhomogeneous plasma. In: Reviews 

of Plasma Physics, Vol . Ill, edited by M. A. Leontovich, Consultants 
Bureau, N. Y. , 1967. 

Mikhailovskii, A. B., and A. M. Fridman, Drift waves in a finite pressure 
plasma, Sov. Phys. JETP , 24 , 965, 1967. 

Parks, G. K. , F. V. Coroniti, R. L. McPherron, and K. A. Anderson, Studies 

of the magnetospheric substorm, 1, Characteristics of modulated electron 
precipitation pccurripg during auroral substorms, J . Geophys . Res , , 


1685, 1968a. 



-42- 


Parks, G. K., T. W. Lezniak, and J. R. Winckler, The Magnetospheric Substorm: 

A fundamental mode of production of energetic electrons in the outer 
Van Allen Belt, paper 6-10. International Symposium on the Physics 
of the Magnetosphere, Washington D. C. , Sept. 1968. 

Petschek, H. E., and C. F. Kennel, Tail flow, auroral precipitation and ring 
currents, Trans. Amer, Geophys. Union , 47 , 137, 1966. 

Rosenbluth, M. N., N. A. Krall, and N. Rostoker, Finite Larmor radius 

stabilization of weakly unstable confined plasmas. N uclear Fusion Suppl. , 
pt. 1, 143, 196$!. 

Rosenbluth, M.N. , Microinstabilities. In: Plasma Physics, International 

Atomic Energy Agency, 1965, p. 485. 

Rosenbluth, M. N. , Low frequency limit of interchange instability, Phys, 
Fluids , _n, 869, 1968. 

Rukhadze, A. A., and V. P. Silin, Method of geometrical optics in the 

electrodynamics of an inhomogeneous plasma, Sov. Phys. Uspekhi , 7, 209, 
1967. 

Rutherford, P. H. , and g. A. Frieman, Drift instabilities in general 
magnetic field configurations, Phys. Fluids , 11 , 569, 1968. 

Stix, T. H., The Theory of Plasma Waves. McGraw-Hill, N.Y., 1962. 

Stringer, T. E., Low frequency waves in an unbounded plasma. Plasma Phys. 

(J. Nucl. Energy, Part C) , 5, 89, 1963. 

Swift, D. Vi. , Anew interpretation of long period micropulsations, 

J. Geophys. Res. , 72 , 4885, 1967. 

Vasyliunas, V. M., A survey of low energy electrons in the evening sector 
of the magnetosphere with OGQ 1 and 0G0 3, J. Geophys. Res. , 73 , 

2839, 1968a. 



-43- 


Vasyliunas, V. M. , Low energy electrons on the day side of the magnetosphere, 
J. Geophys, Res. , 73 , 7519, 1968b. 

Wilson, C. R. , Conjugate three dimensional polarization of high latitude 
micropulsations from a hydromagnetic wave-ionospheric current model, 

J, Geophys. Res. , 71 , 3233, 1966. 



